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ABSTRACT 

Molecular line emission from protoplanetary disks is a powerful tool to constrain their physical and chemical 
structure. Nevertheless, only a few molecules have been detected in disks so far. We take advantage of the 
enhanced capabilities of the IRAM 30m telescope by using the new broad band correlator (FTS) to search for 
so far undetected molecules in the protoplanetary disks surrounding the TTauri stars DM Tau, GO Tau, LkCal5 
and the Herbig Ae star MWC 480. We report the first detection of HC 3 N at 5cr in the GO Tau and MWC 480 
disks with the IRAM 30-m, and in the LkCa 15 disk (5 <x), using the IRAM array, with derived column densities 
of the order of 10 12 cirT 2 . We also obtain stringent upper limits on CCS (N < 1.5 x 10 12 cirT 3 ). We discuss 
the observational results by comparing them to column densities derived from existing chemical disk models 
(computed using the chemical code Nautilus) and based on previous nitrogen and sulfur-bearing molecule 
observations. The observed column densities of HC3N are typically two orders of magnitude lower than the 
existing predictions and appear to be lower in the presence of strong UV flux, suggesting that the molecular 
chemistry is sensitive to the UV penetration through the disk. The CCS upper limits reinforce our model with 
low elemental abundance of sulfur derived from other sulfur-bearing molecules (CS, H2S and SO). 
Keywords: circumstellar matter — protoplanetary disks — stars: individual (DM Tau, LkCa 15, GO Tau, MWC 
480) — Radio lines: stars 
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Understanding the formation of planetary systems requires 
an in-depth study of the initial conditions, i.e. the structures of 
the protopl anetary disks. Sev eral theoretical works have been 
done (e.g. Bergin et al. 2007, and references therein), leading 
to the c urrent picture of a flared disk consisting of th ree lay- 
ers (e.g. Ivan Zadelhoff et alJl200lT : iBergin et al.ll2007l) . How- 
ever, key model parameters such as gas density and tempera- 
ture remain so far poorly constrained by observations due to 
the limited sensitivity and spatial resolution of current instru- 
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merits. Observations of molecular lines have proven to be an 
excellent tool to study the physical and che mical structure and 
the dynamics of pro t oplanetary disks (e.g.lDutrey et al.lll997b 
Kastner etall 119971; IGuilloteau & Dutreyl 11998b IPietu et al] 
20071: lOi et alJ 120081: iSemenov & Wiebd HoTl lOberg et alJ 

Depending on the molecule and transition, the observa- 
tions trace different physical conditions and, therefore, sample 
chemically different regions in the disks. In gas-rich proto- 
planetary disks the C + emission likely comes from the ion- 
ized upper part of dis ks ( Seme nov et alj|2004t iJonkheid et al. 
2007 1 : [Chapillon et all 1201 Ot iPanic et al.l 120 10t iKamp et al. 
201 It Bruderei et al.l 120121) . Atomic species such as Ol 
have also been detected by the Herschel satellite but its ori- 
gin (warm atmosphe re, disk wind or jet?) is still debated 
(Math ews et alJ2 010). CO and its main isotopomers are char- 
acteristics of the molecular layers which are located up to 
a few scale heights only above the mid-plane (Darto lset alJ 
120031) . Near the disk plane, the gas is very cold, so that 
molecules are expected to be deplete d on grains. Such re- 
gions may be only traced by H2D + (IChapillon et alJ 120121 : 
Oberg et al. 201 lb). To retrieve useful information and to re- 
construct the 3D (physical and chemical) structure of disks, 
molecular observations need to be compare d with chemi- 
cal models dedicated to p rotoplanetary disks (ISemenov et alJ 
l2010t IVasvunin et al.ll201 ll) . Such studies help characterise 
the dominant processes in the disk leading to planet forma- 
tion (e.g., grain growth, molecule formation and destruction, 
gas-to-dust ratio evolution). 

So far, in the mm/submm range which characterizes the 
bulk of the gas disk, only a few molecules have been firmly 
detected in protoplanetary disks: CO and its main iso- 
topomers ( 13 CO and C ls O), HCO + and H 13 CO + , DCO + , 
H?CO, H 7 0, CS, C?H, N 7 H + , HCN, HNC, CN and DCN 
dDutrev et all 119971: Ivan Dishoeck et all 120031: IDutrey et all 
120071: lOi et al. 2008; Hen ning et ai1l2010HOberg et alj|201 lab 
Hogerheij de et al.ll201 lb lOberg et alJl2012l) . These are sim- 
ple, light (mass number < 44) molecules. Heavier or more 
complex molecules remain undetected due to their low abun- 
dances and the lack of sensitivity of current instruments. 

Within the framework of the CID ("Chemistry In Disks") 
collaboration, we take here advantage of the improved perfor- 
mance of the IRAM-30m telescope to carry out a molecular 
survey of four well known large (R ol „ > 500 AU) protoplan- 
etary disks. The main driver of our observational project was 
to search for heavier molecules. We report in this observa- 
tional paper the detection of the new molecule HC3N, in the 
three disks surrounding the T Tauri stars GO Tau and LkCal5 
and the Herbig Ae star MWC 480. We also present the up- 
per limits obtained on CCS after a deep integration. A more 
complete chemical analysis will be presented in a forthcom- 
ing paper (Wakelam et al. in prep.). 

In Section 2, we present the observations and the methods 
used for the data reduction. Results (derivation of the best-fit 
model and column densities) are described in Section 3. In 
Section 4 the implications of our new results are discussed. 
Summary and conclusion are presented in Section 5. 

2. SOURCE SELECTION AND OBSERVATIONS 

We focus our survey on well-known molecular-rich 
protoplanetary disks: MWC 480, LkCal5 and DM Tau 
(IDutrey et alJ 1997b IGuilloteau & Dutreyl 1998blHenning et al] 
120 lOt lOberg etalj|20 10ft . and GO Tau dSchaefer et al]|2009ft . 



MWC 480 is an He rbig Ae star with a stellar mass of 1.8 M 
dSimon et al.ll2000l) . while DM Tau and LkCal5 are TTauri 
stars of mass 0.5 and 1 .0 M n , respective ly. Their disk struc- 
tures are well known (IPietu et al.l 120071) . A large CO disk 
(~ 900 AU) around th e 0.5 M Q star GO Tau was reported by 
Schaefer et al. (2009). This di sk is similar in ma ny aspects 
to that surrounding DM Tau dDutrev et alJ 1201 ll) . Sulfur- 
bearing molecules have been recen tly searched for in these 
four disks by the CID collaboration dDutrev et al.ll201 ll here- 
after CIDV). 

Observations were carried out in August 2011 with the 
IRAM 30-m telescope under quite good weather conditions 
using the wobbl er switch and , follow ing our standard observ- 
ing mode dDutrev et al]| 19971 I20TT1) . We used the new 30-m 
heterodyne receiver EMIR. All the lines were simultaneously 
observed using one frequency setup thanks to the capability of 
the new broad-band spectrometer (FTS). The setup covered 
the HC 3 N J=10-9, J=12-ll and J=16-15 rotational lines (at 
frequencies 90.978, 109.173 and 145.560 GHz) and the CCS 
(7,7-6,8), (7,8-6,7) and (11,12-10,11) transitions at 90.686, 
93.870 and 144.244GHz respectively. The frequency reso- 
lution was 200 kHz (i.e. 0.6 - 0.4 km.s -1 ). Frequent pointing 
and focus checks were done on Jupiter, Mars and the quasar 
0528+134. IK Tau/NML Tau served as line calibrator. We 
spent 5 to 12 hours on each source and the T sys was 1 10-130 
K at 3 mm and 135-200 K at 2 mm. The observations were 
converted in unit of flux density using the IRAM recommen- 
dation^. 

We also obtained low spectral resolution (8 km.s -1 ) images 
of LkCal5 and GO Tau with 3.5" angular resolution using the 
IRAM Plateau de Bure interferometer, covering the frequency 
band 144.65 - 148.35 GHz. Observations were performed on 
30- Jul and 14-Sep 2011 in the CD configuration. Standard 
phase and amplitude calibration were applied, and the flux 
calibration used MWC 349 as a reference. 

The data was reduced using the GILDAS0 package. 

3. RESULTS AND COLUMN DENSITY DERIVATION 

Figure Q] is a montage showing the 30-m observational re- 
sults. We emphasize the sensitivity by summing up all ob- 
served transitions of each molecule. The averaged spectra are 
presented in the bottom panels. CCS is not detected but we get 
good upper limits (see TableQ]). The first cyanopolyyne HC3N 
is convincingly detected in two disks at the level of 5<x (GO 
Tau and MWC 480) while it is marginal in LkCal5 (4 <x) on 
the integrated spectra and undetected in DM Tau. MWC 480, 
DM Tau and LkCal5 are located in "holes" of their parent 
molecular clouds in regions devoid of CO. This is not the 
case of GO Tau; ISchaefer et ail (|2009) reported from PdBI 
data that the red part of the CO J=l-0 and J=2-l lines are hid- 
den by the presence of CO at the same velocity. We carefully 
checked if confusion can affect th e HC3N flux density in GO 
Tau. From Schaefer et al.l d2009l) and 13 CO data obtained si- 
multaneously with the HC3N observations, we estimate which 
velocity channels could suffer from contamination. While the 
13 CO data still present some confusion with the CO ambient 
cloud, this is not significant for the HC3N lines. Confusion, if 
any, would rather absorb the red-shifted part of the spectrum 
than lead to excess emission, because the critical densities for 
substantial excitation of the J=ll to J=16 rotational levels is 
high. Any analysis using only the unconfused half of the disk 

1 http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 

2 http://www.iram.fr/IRAMFR/GILDAS/ 
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Figure 1. Observations of HC3N and CCS in the four disks. Top panel: HC3N. Bottom panel: CCS. In both panels, the last line is the average of the three 
observed transitions. The best fits or the 3<x upper limits of the best model are also superimposed on the individual lines. 
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Figure 2. Spectrum of HC3N and H 2 CO towards LkCal5 obtained at 3.5" 
resolution with the IRAM array. The two detected lines are indicated. 

velocity profile would lead to an increase of the HC3N col- 
umn density estimate. For GO Tau, only an upper limit of 
< 0.1 Jy.km/s can be reported for H3CN (16-15) from the 
PdBI observations, compatible with the upper limit given in 
Table Q] 

For LkCa 15, the blue-shifted part appears somewhat 
stronger than the red-shifted one. Spectra of 13 CO taken si- 
multaneously do not show such a behavior, so this cannot be 
ascribed to mis-pointing. Note that this effect is however at 
the Signal to Noise limit. Figure [2] shows the spectrum to- 
wards LkCal5 at 3.5" spatial resolution obtained with the 
PdBI, after removal of the 20 mJy/beam continuum emission 
given by the wideband correlator. The velocity scale refers 
to the rest frequency of the J= 16-15 transition of HC3N at 
145.560938 GHz. The line is clearly detected (S/N > 5), 
with a total line flux around 0.15 Jy.km/s, and is about a fac- 
tor 1.5 weaker than the 2o,2 - lo,i transition of para H2CO at 
145.602953 GHz also visible in this spectrum. These obser- 
vations clearly show that the emission detected with the 30-m 
is indeed coming from the circumstellar disk. 

We retrieve the molecular column densities from the data 
using DiskFit (a de dicated radiative transfer code for proto- 
planetary disks, see iPietu et al.ll2.Q07h and following the stan- 
dard method described in detail in CIDV where the physi- 
cal disks models are parameterized by power laws. The disk 
geometry and kinematics are consistently incorporated in the 
model. In particular, we properly take into account the ve- 
locity gradient due to Keplerian rotation. At a given veloc- 
ity, only one well defined portion of the disk emits and hence 
contributes to the observed column density. All parameters 
were fixed (see Table |2]i except the molecular surface densi- 
ties (see Section 2.2 of CIDV for more information). Since 
we have observed three lines of the same molecule, both for 
CCS and HC3N, the best disk model is obtained b y fitting si- 
multan eously these three lines, as we proceeded in lPietu et al] 
(2007) for CO. Note that the best model is obtained using 
a natural weighting (system temperature only) of each line. 
As in the analysis of the sulfur data (CIDV), we assume for 
the surface density power law versus radius with an exponent 
p = 1.5 (I(r) = 2300 x (r/300AU) -1,s and for the tempera- 
ture T(r) = rujoWlOOAlTr 9 with q = 0.4. As a radius of 
reference, we use 300 AU because it is representative of the 
area where the unresolved and beam diluted emission arises. 
A smaller value would be unrealistic because of the beam di- 
lution effect, a larger value would correspond to a disk area 
where lines would be too optically thin. For GO Tau, the 
source with the highest S/N ratio, we also tried to fit simulta- 



neously p with E300 and found p = 1.4 + 0.3. Table|2] adapted 
from CIDV, summarizes the physical disk properties used to 
retrieve the best fit disk model. The derived surface densities 
at 300 AU are presented in Table[3] The line profiles obtained 
with the best disk models are overplotted in Fig[T]on the ob- 
servations. The best fits do not follow the individual spectra 
because it is the result of a simultaneous fit of the three lines 
giving more weight to the best S/N ratio spectrum. 

4. DISCUSSION 
4. 1. Chemical models of protoplanetary disks 

Chemical modeling is essential to retrieve the vertical 
structure of the disk from the observable column densities. 
Generic chem ical models of cir c umste lla r disks have been 
publis hed bv lAikawa & Herbstl (1 19991) . IWillacv & Langerj 
(2000), lAikawa et al l (120021) Ivan Zadelhoff et al l 
d2003l) lAikawa & Nomural (120061). IWillacv et al l (120061) 
and lFogeletal.1 (120111) ; ISemenov & Wiebd (120111) ; 
IVasvunin et al.l (1201 llK Each of these works make dif- 
ferent assumption and hypothesis. lAikawa & Herbstl (119991) 
consider gas phase chemistry, along with sticking and 
desorption, and assume a vertically isothermal disk model 
following the "mini mum mass solar nebula" extended to 800 
AU. The study from Willacy & Langer (2000) use the th ermal 
structure derived from the IChiang & Goldreichl ( 1 19971) two- 
layer approxima tion, as well as enh a nced photode sorption 
yields following IWestlev et al l (119951) . lAikawa et ail (120021) 
use a D'Alessio disk model including vertical temperature 
gradients and self-consistently variable flaring, although 
dust and gas are assumed to be fully thermally coupled. 
All three models use a 1 + 1D approximation for radiative 
transfer, in which the stellar UV is attenuated along the line 
of sight to the star, while the ISRF impacts is otropically 
on the disk surface. Ivan Zadelhoff et a l. (2003) improved 
the UV treatment by using a 2-D radiative tran sfer code to 
solve for the UV field inside the disk, as did Fogel et al. 
d201 ll) . Shieldi ng by H? is treated in a n approximate 
way, however. lAikawa & Nomural (120061) expanded the 
models furth e r by c on sidering the effect o f grai n growth as 
Woitk e"eTaT] (120091) . IWoitke et al.l (120091) and [ Fogel et al l 
(120111) investigated the effec ts of dust settlin g, Fogel et all 
mUb . IWillacv et all (120061) and iWalsh et al.l (120101 120121) 
also include grain-surface chemistry, which may be important 
in such environ ments. ID tur bulent diffusion is included in 
the model from IWillacv et al.l (2006), and 2D in the models 
from ISemenov et all (120061) and ISemenov & Wiebd (120111) . 
IVasvunin et al.l d2011l) have investigated dust coagulation, 
fragmentation, sedimentation, and turbulent stirring.. 

Making comparisons with our observations is difficult for 
several reasons. All these chemical models are not necessar- 
ily tailored to the physical conditions relevant to the specific 
sources we have observed. The chemical networks used are 
different from one code to another and the large uncertainties 
on reaction rates ma y explain some intrinsic differences (e.g. 
iDaranlot et al.l l2012). In many cases the initial physical and 
chemical conditions are also different. 

4.2. The Nautilus models 

Following our previous studies, we performed a chemical 
modeling using Nautilus, a gas-grain chemis try model 
adapted for the disk physics (Her sant et al.| [2009). Nautilus 
computes the abundances of 460 gas-phase and 195 surface 
species as a function of time using the rate equation method 
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Table 1 

IRAM 30-m integrated line fluxes for HC3N lines and 3<x upper limits for CCS (Jy.km.s -1 ) 



Sources 


LkCa 15 


MWC480 


DMTau 


GOTau 


HC3N (16-15) 


0.10 ±0.03 


0.090 ± 0.029 


< 0.06 


<0.08 


HC 3 N (12-11) 


< 0.09 (0.034 ±0.031) 


0.132 ±0.040 


< 0.09 


0.102 ±0.024 


HC 3 N(10-9) 


< 0.06 (0.025 ± 0.022) 


0.060 ± 0.026 


< 0.06 


0.067 ± 0.018 


CCS (11,12-10.11) 


< 0.36 


< 0.11 


< 0.08 


< 0.09 


CCS (7,8-6,7) 


< 0.05 


< 0.08 


< 0.05 


< 0.05 


CCS (7,7-6,8) 


< 0.06 


< 0.08 


< 0.06 


< 0.05 



Note. — Integrated fluxes are derived from gaussian fits with fixed central velocities (V sys given in Tab|5J and line-widths corresponding to Keplerian motions 
with AV = 1.4 km.s-' for GO Tau and DM Tau, AV = 4.1 km.s" 1 for MWC480 and AV = 3.3 km.s" 1 for LkCal5. 



Table 2 

Disk physical parameters used to derive the best fit models 



Source 


LkCa 15 


GOTau 


DM Tau 


MWC 480 


inclination (") 


52 


51 


-32 


38 


P.A.(°) 


150 


112 


65 


57 


Vj VIf (km.s _1 ) 


6.30 


4.89 


6.04 


5.10 


Vioo(km.s _1 ) 


3.00 


2.05 


2.16 


4.03 


^.(km.s" 1 ) 


0.2 


0.2 


0.2 


0.2 


Tioo(K) 


15 


15 


15 


30 


R,„,(AU) 


45 


1 


1 


1 


R ou ,(AU) 


550 


900 


800 


500 



Note. — Following Dutrey et al. 12011). RA. is the position angle of the 
disk rotation axis, i is the inclination, V syst is the systemic velocity and <S„ the 
turbulent line-width component. The velocity laws (V 100) are Keplerian (as 
shown by Pietu et al. 12007) and Schaefer et al. 1 2009)) 



(Hasegawa et alj [19921) and the chemical network contains 
4406 gas-phase reactions and 1733 reactions involving 
grains, including adsorption and desorption processes and 
grain-surface reactions (we assume 4580 K for the energy 
of desorption of CCS and HC3N). The gas-phase network is 
regularly updated (both by optimizing the reaction rates and 
by adding new reactions) according to the recommendations 
from the KIDA0 experts. In CIDV, we used Nautilus to per- 
form a full chemistry modeling and only discussed the sulfur 
chemistry. The initial abundances were obtained computing 
the chemical composition of the parent cloud as discussed 
in CIDV. We compare here our new observational results 
with our best model from CIDV (model C). For the four 
disks, the physical conditions and elemental abundances used 
are summarized in Fig. 2 and Table 4 of CIDV, respectively. 
Model C was obtaine d assuming a C/O ratio of 1 .2 (following 
iHincelin et al.ll20T11) . a sulfur and nitrogen abundance with 
respect to H of 8 x 10~ 9 and 6.2 x 10~ 5 dJenkinsi [20091) . 
respectively and a grain size of 0. 1 /mi with an initial cloud 
density of 2 x 10 5 H.crrT 3 and a cloud age of 1 Myr. Our 
results are shown in Figj3] At 1 Myr, we have an abundance 
of 10~ 10 for HC 3 N and 4 x 10~ 13 for CCS in the gas-phase 
and of 2.5 x 10~ 12 for HC3N and 4 x 10" for CCS, on 
grains. The main conclusion from CIDV was that although 
our chemical model reproduces the SO and CS column 
densities reasonably well, it fails to reproduce the upper 
limits obtained on H2S by at least one order of magnitude. 
This suggests that a fraction of Sulfur may be depleted in 
mantles or refractory grains. At the high densities and low 
temperatures encountered around disk mid-planes, H2S may 
remain locked onto the grain surfaces and react to form other 
species preventing desorption of H2S. 



4.3. Comparison with observations 
3 http://kida.obs.u-bordeauxl.fr 



In protoplanetary disks, only four N bearing molecules have 
been previously detected. NH3 is not detected but N2I-U 
is observed with an abundance relativ e to H 2 of ~ 10" 12 
dDutrev et al. | |2007b lOberg et al.ll20lTah . CN is easy to detect 
( Dutre v et al.ll 19971) and likely results from a chemistry which 
is more sophisticated than the simple photo -dissociation of 
HCN under the stellar and ambient UV fields (Chapillon et al. 
120121) . As an example, in the DM Tau disk, the observed sur- 
face densities of CN and HCN at 300 AU are 3.5 • 10 13 cirT 2 
and 6.5 • 10 12 cm" 2 , respectively . HNC has been reported 
in DM Tau only by iDutrev et all (fl997l) . HC 3 N is the fifth 
nitrogen-bearing molecule detected in protoplanetary disks. 
We observe this molecule in MWC 480, GO Tau and LkCa 
15. Figj3] presents the modeled and observed surface densi- 
ties of HC3N at 300 AU and Fig|4]presents the modeled abun- 
dances. DM Tau and GO Tau are assumed to share the same 
disk physical structure and therefore the chemical predictions 
are identical. 

Figf5] (left panel) shows that our non detection of CCS 
remains, at the 3<x level, in very good agreement with the 
model C from CIDV. Furthermore, our CCS upper limits are 
incompatible with the other models we explored in CIDV. 
By increasing the elemental abundance of sulfur by a fac- 
tor 10 (models A and B), we would produce ten times more 
CCS with respect to H2. Figf3] (right panel) also shows that 
CCS is expected to be formed around 1.5 scale-heights in 
the outer disk (R~ 300 AU ) of both TTauri and Herbig 
Ae stars. Sem enov & Wiebel (1201 11) have calculated typical 
column densities of gaseous CCS at 300 AU to be about 
10 11 cirT 2 and 710 11 cirT 2 in the laminar and 2D-mixing 
models of the DM Tau disk. The CCS layer in these mod- 
els is also located around 1 pressure scale height from the 
midplane. 

Contrary to the observed behavior of sulfur-bearing 
molecules (CS and H2S), whose surface densities are of the 
same order of magnitude in both sources, HC3N is only de- 
tected in the disk surrounding GO Tau. In DM Tau, the 3cr 
upper limit is a factor 4 lower than the detection in GO Tau. 
Some intrinsic differences in the physical/chemical properties 
have to be explored to understand the origin of such a differ- 
ence. 

4.4. Results 

For the TTauri disks, our best chemical model predicts 
HC3N surface densities which are typically two orders of 
magnitude larger than the observed one. Our chemical model 
show that HC3N is found in a layer around 1-2 scale height 
in the TTauri disks (see Fig|3j. In this layer, the density is of 
the order of ~ 3 x 10 6 crrT 3 , i.e. comparable or sli ghtly above 
the ex pected critical densities of the transitions dWernli et al.l 
12007b . especially for J=16-15 transition. As we assumedLTE, 
this may affect our determination of the surface densities by a 
factor of at most a few, but not as large as ~ 100. 
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Table 3 

Derived and predicted surface densities and ratios 







^300 


(cm 2 ) 














HC 3 N 


CCS 


HC 3 N/HCN 


HC3N/CN 


Souce 


Derived 


Predicted 


Derived 


Predicted 


Derived 


Predicted 


Derived 


Predicted 


LkCa 15 


8 ±2- 10 u 


5.2- 10 lj 


< 1.4- 10 12 


2.9- 10" 


0.075 


0.12 


0.014 


0.011 


GOTau 


13 ± 2- 10" 


4.4- 10 13 


< 1.2- 10 12 


3.7- 10" 




0.08 




0.02 


DMTau 


<3.5- 10" 


4.4- 10 13 


< 1.1- 10 12 


3.7- 10" 


< 0.05 


0.08 


<0.01 


0.02 


MWC480 


6± !• 10" 


6.4- 10" 


<0.9- 10 12 


3.1- 10" 


0.55 


0.05 


0.057 


0.001 



Note. — Col. 2 and 4: Surface densities derived at 300 AU from the 30-m data using Diskfit (see section [2] for details). Col. 3 and 5: 
Predicted surface densities at 300 AU using Nautilus from CIDV best chemical model (see section|4]for details). Note that the model is the 
same for DM Tau and GO Tau. Col. 6-9: Surface densities ratios at 300 AU. "Derived ratios" are calculated using the CN and HCN surface 
densities in LkCa 1 5 , DM Tau and MWC 480 from Chapill oiTeTall i20l2l) . 



The chemical model however appears in good agreement 
for the disk surrounding the Herbig Ae star MWC 480 in 
which the predicted abundance is much smaller than in the 
other sources. We checked that the difference is due to the 
UV flux, significantly stronger in the case of MWC 480, 
which photo-dissociates HC3N. The HC3N abundance 
appears better reproduced in the presence of a strong UV 
field. Enhancing the UV penetration by changing vertically 
the dust grain properties (grain growth and vertical settling) 
may decrease the predicted HC3N column densities even in 
sources with low or medium UV fluxes. The stronger UV 
flux also explains that most of the HC3N is formed at an 
altitude which is lower in the Herbig Ae disk (Z/H ^ 1) than 
in TTauri disks (Z/H ^ 1.5, see Figf3] right panel). As are 
result, the observed HC3N transitions should be more easily 
thermalized in the Herbig Ae disk. 

Investigating the vertical abundance variation is also very 
intersting. Figure [4] shows the predicted vertical abundances 
of HC3N and CCS, in the gas-phase and at the surface of 
the grains at 300 AU in the four disks. These two molecules 
are present in a layer between 1 and 3 scale heights in DM 
Tau and at lower height in the two other disks where the 
stellar UV flux is stronger (410, 2550, 8500 xo for DM 
Tau, LkCa 15 and MWC 480 respectively, see CIDV). In all 
models, CCS is formed in the gas-phase only and can only 
be destroyed on grain surfaces, after being accreted. In DM 
Tau, the surface abundance of CCS is smaller than in LkCa 
15 because more complex S-bearing species such as C3S 
are formed. In LkCa 15, the formation of such molecules 
is prevented by the stronger UV field. In MWC 480, the 
CCS abundance in the ices is even smaller than 10~ 12 . The 
gas phase abundance of HC3N for DM Tau is similar to the 
one in LkCa 15, whereas the abundance in the grain surface 
is much larger. The smaller surface abundance of HC3N in 
LkCa 15 is due to direct photo-dissociations by UV photons. 
The situation in MWC 480 is very different, with a very 
large abundance of HC3N in the ices and a very small one in 
the gas. Although the UV field is stronger than in the two 
other sources, the dust temperature is larger and allows the 
diffusion of larger radicals on the surfaces. This enhanced 
diffusion produces more HC3N on grains. 

In Table [3] we co mpare our results with those from 
Chapi llon etatl (12012b . If our prediction of the HC3N col- 
umn density in MWC 480 is in good agreement with the ob- 
servation, it fails to reproduce the HC3N/CN and HC3N/HCN 
ratios. For the T Tauri disks, the ratios are in roughly good 
agreement with the observations while the predicted absolute 
values are off by two orders of magnitude. This is not surpris- 



ing since this model was aimed at reproducing sulfur-bearing 
molecule chemistry. Using a lower elemental abundance 
of nitrogen would produ ce less N bearing molecules, (see 
ISemenov & Wiebal201 lh . Moreover, some neutral-neutral re- 
actions involving atomic nitrogen are currently being exper- 
imentally studied dBergeat etal.ll2009t iDaranlot et al] 120 1 lL 
120121) . As soon as these new experimental rates are known, 
a more det ailed study of nitroge n-chemistry in disks will be 
performed ( Wakela m et al.ll2012l) . We also note that the ob- 
served CN/HCN ratio is of the order of 5 in the T Tauri stars 
disks, and of 10 in the disks of Herbig stars. 

The observed HC3N/HCN ratio in the three disks is 
very similar to that observe d in the cold molec ular core 
L134N (HC3N/HCN ~ 0.07, lDickenslFaT]|2000T) and stiU 
reasonably close to the cometary value within a factor ~ 2 
(in H ale-Bopp, HC3N is observed with a ratio HC3N/HCN 
~ O.l.lBockelee-Mo rvan et alJl2000h . However, the observed 
HC3N/CN ratio of both Herbig Ae and TTauri disks is 
far from the value observed in L134N (HC3N/CN ~ 1, 
iDickens et al.l 120 00), suggesting that the later ratio is more 
affected by photo-dissociation. 

These observations can be used to place the detectability of 
complex molecules in perspective with ALMA. Our typical 
3cr limit is of order 50 mJy.km.s~ 1 for a single line (Table[T|), 
but we gain V3 by the multiplex factor due to wide frequency 
coverage. By comparison, ALMA reaches 1 mJy in 1 hour of 
integration time for 1 km.s 1 resolution at 90 GHz, i.e. \<j flux 
sensitivity of 1.4 mJy.km.s~ 1 . ALMA should greatly improve 
the detectability of such molecules, but imaging at high (< 
0.5") resolution would require many (more than 10) hours. 

5. SUMMARY 

Using the 30-m IRAM telescope, we report the first detec- 
tion of HC3N in three disks surrounding young low-mass stars 
and present deep upper limits on the CCS molecule. Our re- 
sults allow us to conclude: 

• The CCS upper limit confirms our previous modeling of 
sulfur-bearing molecule which favors a relatively low 
elemental abundance of sulfur (8 • 10~ 9 /H). 

• The observed column densities of HC3N are typically 
two orders to magnitude lower than the predictions. 

• A strong UV flux decreases the abundance of HC3N 
suggesting that a better UV penetration either due to 
significant grain growth or dust settling may signifi- 
cantly affect the HC3N surface density, even when a 
moderate UV flux impacts the disk. 
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Figure 3. Predicted column densities for HCjNand CCS using Nautilus. 
Left panel: observed and modeled surface densities. Right panel: Cumu- 
lative column densities (from mid-plane to atmosphere) in function of Z/H 
(left axis) and Z (right axis) at 300 AU. Top panel: MWC 480, medium panel: 
LkCal5, Bottom panel: GO Tau and DM Tau. The (cumulative) column den- 
sities (right panel), whose asymptotic values are half of the surface densities 
(left panel). 

• More generally, our observations confirms that complex 
molecules remain difficult to detect in disks. 
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